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1 INTRODUCTION 



ABSTRACT 

We present a new unbiased search and analysis of all B stars in the solar neighbourhood 
(within a volume of 400 pc diameter) using the Arivel data base to track down the 
remains of the OB associations, which hosted the supernovae responsible for the Local 
Bubble in the interstellar gas. We find after careful dereddening and by comparison 
with theoretical isochrones, that besides the Upper Scorpius the Upper Centaurus 
Lupus and Lower Centaurus Crux subgroups are the youngest stellar associations in 
the solar neighbourhood with ages of 20 to 30 Myr, in agreement with previous work. 
In search for the "smoking gun" of the origin of the Local Bubble, we have traced the 
paths of the associations back into the past and found that they entered the present 
bubble region 10 to 15 Myr ago. We argue that the Local Bubble began to form then 
and estimate that 14 to 20 supernovae have gone off since. It is shown that the implied 
energy input is sufficient to excavate a bubble of the presently observed size. 

Key words: solar neighbourhood - open clusters and associations: individual: 
Sco 0B2 - ISM: individual: local bubble - Local ISM 
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The Local Bubble (LB), a low-density X-ray emitting cav- 
ity deficient of Hi, is our Galactic habitat. Yet, until re- 
cently, its origin remained mysterious. It was conjectured 
to be the result of one or several supernova explosions (e.g. 
Cox & Anderson 1982, Innes & Hartquist 1984, Smith & 
Cox 2001), but firm evidence was lacking, as no OB asso- 
ciation was found within its boundaries, extending about 
200 pc in the Galactic plane, and 600 pc perpendicular to 
it, but inclined by about 20° with respect to the axis of 
Galactic rotation, similar to Gould's Belt (cf. Lallement et 
al. 2003). Further problems arose, since the X-ray and EUV 
spectra measured in the Wisconsin Survey, by ROSAT PSPC, 
DXS, XQC and EUVE were severely at odds with a ther- 
mal hot plasma in coUisional ionization equilibrium (CIE) 
as was pointed out by Jelinsky et al. (1995), Sanders et 
al. (2001), McCammon et al. (2002). Most recently Hurwitz 
et al. (2005) analyzed CHIPS data and found an extremely 
low emissivity of EUV iron lines. The underabundance of 
soft X-ray lines can be naturally explained if the plasma is in 
a state of delayed recombination (Breitschwerdt & Schmut- 
zler 1994, Breitschwerdt 2001), but a high resolution nu- 
merical hydrodynamical evolution model is needed to bet- 
ter constrain non-equilibrium models. Spectral discrepancies 
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between models and observations can be alleviated if there 
is a substantial contribution from very local sources, such as 
the Earth's exosphere (Freyberg 1998) or charge exchange 
reactions between solar wind ions (SWCE) and heliospheric 
gas (LaUement 2004). At present it is unclear what fraction 
can be attributed to these very nearby sources, although 
there is fairly robust evidence that even in the extreme case 
of all of the X-ray emission being due to SWCE in a cer- 
tain direction, a substantial LB fraction remains, especially 
perpendicular to the disc. For further details on LB proper- 
ties we refer to the review of Breitschwerdt (2001) and the 
conference proceedings The Local Bubble and Beyond (Bre- 
itschwerdt et al. 1998). 

All these shortcomings have led several authors to spec- 
ulate that if the LB is not a classical superbubble, but rather 
an appendix of the neighbouring Loop I superbubble, which 
was expanding into an interarm region between the Sagit- 
tarius and the Perseus spiral arms of the Galaxy (Bochkarev 
1987, Frisch 1995). However, the existence of a "wall" be- 
tween the two bubbles, showing up in absorption of soft 
X-rays in ROSAT PSPC images (Egger & Aschenbach 1995) 
renders this scenario not very plausible. 

The search for the "smoking gun" of the origin of the 
LB proved partially successful by discovering that moving 
groups of young stars in the solar neighbourhood could pro- 
vide an adequate number of SN explosions while crossing the 
path of the LB. Berghofer & Breitschwerdt (2002, henceforth 
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BB02) calculated the trajectory of the Pleiades subgroup Bl 
backwards in time, and found that 19 SNc could have ex- 
ploded between 10-20 Myr ago in the region that is occupied 
by the LB. The remaining stars of Bl are now part of the 
Scorpius Ccntaurus OB association. It could be shown that 
this is in good agreement with the size of the LB and the 
present soft X-ray emissivity. A similar analysis was carried 
out by Maiz-Apellaniz (2001), who calculated backwards in 
time the trajectories of Sco Cen subgroups and claimed that 
about 6 SNo that went off in the Lower Ccntaurus Crux sub- 
group of the Sco 0B2 association 7 to 9 Myr ago formed the 
LB. 

While these analyses represent a major step towards 
the understanding of the origin of the LB they arc not free 
from bias, in particular the assumption that certain stel- 
lar groups should be responsible for the sought SN explo- 
sions. The purpose of this paper is to scrutinize all stars 
that are within a volume of about 400 pc in diameter cen- 
tered around the Sun, and to perform a selection according 
to spectral and kinematical properties. The latter is based 
on three dimensional space velocities of the stars. Thus our 
approach is complementary to studies like by de Zeeuw et 
al. (1999) which are based on proper motions alone. Sartori 
et al. (2003) do include radial velocities when analyzing the 
subgroups of the Sco 0B2 association, but work from a list 
of stars preselected by de Zeeuw et al. (1999). From their 
position in the HR diagram and the turnoff point from the 
ZAMS, we can reliably determine the age of the stars and 
estimate the number of SNe within a defined region, such as 
the LB. 

The paper is organized as follows. In the next section 
we describe our search strategy for the remnants of the OB 

association responsible for the origin of the LB. In section 3 
we discuss the consistency of our findings with the properties 
of the LB as observed today, and present a high-resolution 
3D hydrodynamical simulation of the formation of the LB in 
the local interstellar gas. In the final section we summarize 
our conclusions. 



2 SEARCH FOR NEARBY OB ASSOCIATIONS 

Even though the Sco 0B2 association has been claimed with 

good reason to be responsible for the origin of the LB 
(Maiz-Apellaniz 2001, BB02), a fresh, unprejudiced search 
for the OB association, that might have triggered the for- 
mation of the LB, will improve and harmonize previous 
studies. As starting point we used the Hipparcos catalogue 
(ESA 1997) from which we selected all stars bluer than 
[B — V) < —0.05 with parallaxes larger than 5 milliarc- 
seconds, because otherwise the distances would be too in- 
accurate. Drimmel et al. (2000) find by a comparison with 
the Tycho catalogue that the Hipparcos catalogue is about 
97 percent complete down to V = 7.5 which corresponds at 
a distance of l/(5mas) — 200 pc to an absolute magnitude 
of Mv = 1. Allowing for an extinction of Ay ^ 0.2 we have 
My < 1.2, and we have chosen the colour cut in (B — V) ac- 
cordingly (cf. Figs. 4 and 5). We have obtained this way an 
unbiased complete sample of 762 B stars within a distance 
of 200 pc from the Sun. One of the bluest stars is the BO. 5 
star aCru. The stars of the sample show at this stage al- 
ready a well defined main sequence in the colour-magnitude 



Figure 1. Positions of 610 stars drawn from the Hipparcos cat- 
alogue. The selected stars have colours B — V < —0.05 and for 
each star its radial velocity is known. The X-axis points towards 
the Galactic centre, Y into the direction of Galactic rotation, and 
Z towards the Galactic north-pole, respectively. 



diagram (cf. Fig. 4) . There are a few stars lying several mag- 
nitudes below the main sequence, which are probably sub- 
dwarfs. We have omitted these, My > 3.0 mag, because we 
are interested in young stars. All stars in our list appear in 
the Arivel data base, which combines Hipparcos parallaxes 
and proper motions with radial velocity data collected from 
the literature (Wielen et al., in preparation). We find that 
for 610 stars radial velocities are available. These are accu- 
rate enough to be used for a kinematical analysis (cf. Fig. 3). 
We have tested if the availability of radial velocities is corre- 
lated with the colours of the stars and find that the 152 out 
of 762 stars for which no radial velocity is known lie above 
{B — V) > —0.1 . We discuss the implications of this effect 
on the determination of the expected number of supernovae 
below. For each star spatial X, Y, and Z coordinates and 
the corresponding U, V, and W velocity components have 
been calculated. The coordinates are centered on the Sun 
with X pointing towards l = 0°,b = 0°,Y towards / = 90°, 
& = 0°, and Z towards b = 90°, respectively. Since OB asso- 
ciatious disperse slowly on time scales of the order of several 
10^ years (Blaauw 1964), we have searched for kinematically 
coherent structures in our sample. For this purpose we have 
traced the positions of the stars backward in time over 3-10^ 
years. Stellar orbits have been calculated using the epicyclic 
equations of motion for the stars in the sample (Lindblad 
1959, Wielen 1982) 

X{t) = X(0)-^(l-cos(/tt)) + ^sin(Kt), 



U{t) 



-2B' 
U{0) cos (Kt) 



-2B 



Y{t) = Y{0) + 2A X{0) 



V{Q) sin(/tt), 
-2B I 



(1) 



+ ^F(0)sin(Kt) + ^[7(0)(l-cos(/tt)), 



Vit) = 

m = 



-2B 

K 

W{Q) 



C/(0) sin [Kt) -h V{Q) cos (/tt), 

sin{vt) + Z{Q) cos {vt), 



W{t) = W{Q)cos{vt)- Z{Q)vsm{ut). 

In Eqns. (1) k denotes the epicyclic frequency, k = 

a/— 4f2o_B. A and B are the Oort constants, and fio is the 
angular frequency of the rotation of the local standard of 
rest around the Galactic centre, Oo = Vhsn/Ro- ^ denotes 
the vertical oscillation frequency which is related to the local 
density po by the Poisson equation as v = V AirGpo , where 
G is the constant of gravitation. For the angular velocity of 
the local standard of rest we have adopted a value of Qo 
= 220 km/s/8kpc. The choice of the Oort constants was 
guided by the consideration that they describe in eqns. (1) 
the smooth Galactic gravitational potential. The latter is 
consistent with an essentially flat shape of the local Galactic 
rotation curve, A = -B = Q.o/2 (Feast & Whitelock 1997). 
This must not be confused with determinations of A and B 
using OB stars as, for instance, in the studies of Torra et 
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Figure 2. Positions of the originally selected stars today (pink) 
and S-IO'^ years ago (blue). The Sun is at rest in the diagrams. 
Stars lying in the windows indicated by dashed lines axe identified 
as putative members of the searched for OB association. 



Figure 4. Colour-magnitude diagram of the final sample (236 
stars). Members of the Upper Centaurus Lupus subgroup of 
Sco OB2 are highlighted in orange, Lower Centaurus Crux in yel- 
low, and Upper Scorpius in grey, respectively. 



Figure 3. Present day velocity distribution of the 302 selected 
stars. A second selection is made of the stars lying in the windows 
indicated by dashed lines. 



al. (2000) or Elias et al. (2006). These reflect peculiarities of 
the orbits of the OB stars in Gould's belt related to the ve- 
locities with which they were born, but not the characteristic 
smooth shape of the Galactic potential. For the local density 
we adopt a value of po = 0.1 A4o/pc^ (Holmberg & Flynn 
2004). These parameter values imply k = 0.039 km/s/pc = 
4-10"®yr~^ and u = 0.074 km/s/pc = 7.5-10"*yr"^ In Fig. 2 
we show the positions of the stars today and S-IO'^ years ago. 
Apparently most stars came from directions —90° < I < 90° 
and stayed close to the Galactic midplane. Most of the 610 
stars do not belong to the OB association, which hosted the 
SNc responsible for the origin of the LB, and have space ve- 
locities different from the velocity of the association. Thus 
they are dispersed away into a wide cloud. However, the 
overdense regions in Fig. 2 indicate that there is a consider- 
able number of stars which stayed together. The larger size 
of the overdense regions in the back projected sample com- 
pared to its size today is obviously due to the observational 
errors. The typical accuracy of Hipparcos proper motions is 
about 1 mas/yr which corresponds at a distance of 100 pc 
to a velocity of 0.5 km/s, whereas the accuracy of the radial 
velocities is several km/s. Taken together with an expansion 
velocity of the order of 10 km/s (Blaauw 1964), this implies 
a spreading of the overdensity, which represents the kine- 
matically homogenous group of stars, to a size of roughly 
500 X 500 pc in X and Y. As can be seen from Fig. 2 there 
is an outer shroud of stars which lies at greater distances 
from the core of the overdensity. These must be stars with 
genuinely different space velocities from the kinematically 
homogenous group of stars. We identify this kinematically 
homogeneous group of stars as an OB association and se- 
lect 302 stars lying in the windows indicated as dashed lines 
in Fig. 2. As expected these stars are more or less closely 
related to the Sco 0B2 association. 

In Fig. 3 we show the present day velocity distribution 
of the 302 selected stars. Since the velocity dispersion of an 
OB association is of the order of 10 km/s (Blaauw 1964) 
or even loss (Kamaya 2004), we make a second selection in- 
dicated by windows drawn as dashed lines in Fig. 3. This 
leaves a sample of 236 stars which we analyze in the fol- 
lowing. Hipparcos numbers of these stars are listed in the 
appendix. 

The final sample is shown as a colour-magnitude dia- 
gram in Fig. 4. For this purpose we have cross-identified the 
sample stars in the Geneva photometry data base (Mermil- 
liod, Hauck & Mermilliod 2000) and replaced the {B - V)t 
colours given in the Hipparcos catalogue by {B — V)j colours, 
because they can be then directly compared with theoreti- 
cal isochrones available in the literature. In the colour range, 
which we consider here, B — V given in the Tycho system 
can not be transformed directly to the Johnson system (ESA 



Figure 5. Dereddened colour-magnitude diagram of the mem- 
bers of the Upper Scorpius (grey), Upper Centaurus Lupus (or- 
ange) and Lower Centaurus Crux (yellow) subgroups. The solid 
lines are theoretical isochrones colour coded according to their 
ages. 

1997). The absolute magnitudes have been determined from 
the visual magnitudes given in the Hipparcos catalogue in 
the Johnson system. 

We have compared our sample with the extensive mem- 
bership list of the Sco OB2 association compiled by de Zeeuw 
et al. (1999) who applied a combination of a modified con- 
vergent point method and the so called spaghetti method 
(Hoogerwerf & Aguilar 1999) to Hipparcos data. Of particu- 
lar interest are the membership lists of the subgroups Upper 
Scorpius (US), Upper Centaurus Lupus (UCL), and Lower 
Centaurus Crux (LCC). With only very few exceptions all 
stars in the membership lists, which fulfill our colour selec- 
tion criterion, appear also in our sample, which gives confi- 
dence in our selection procedure. A few stars from our final 
sample could be identified additionally in the membership 
list of de Geus et al. (1989) as members of the subgroups. 
The 79 stars common to both lists are colour coded in Fig. 4 
and listed separately in the appendix. 



3 RESULTS AND DISCUSSION 

3.1 The search for the "smoking gun" 

The colour-magnitude diagram presented in Fig. 4 shows a 
clearly discernible main sequence, which is particularly well 
delineated by the members of the UCL group. The turn- 
off point at the tip is defined by both the members of the 
UCL and the LCC subgroups. Apparently these are together 
with the US subgroup indeed the youngest OB associations 
in the solar neighbourhood (de Geus et al. 1989, Sartori et 
al. 2003) . In order to determine their age we have compared 
the colour-magnitude diagram with theoretical isochrones 
calculated by Schaller et al. (1992) for solar metallicities. 
Fortunately de Bruijne (1999) and Sartori et al. (2003) have 
determined individually for most members of the US, UCL, 
and LCC subgroups, respectively, the extinction and colour 
excess by comparing the observed {V — I)c colours with 
the intrinsic colours of stars of the same spectral type and 
luminosity class. Dereddened data of the 79 stars are shown 
together with isochrones in Fig. 5. We conclude from Fig. 5 
that the ages of the UCL and LCC subgroups lie in the 
range of 20 to 30 Myr, whereas we cannot date the age of 
the US subgroup on the basis of our data. We note that 
this estimate of the ages of the subgroups is nearly twice 
of that of de Geus et al. (1989), who determined an age 
of 11 - 12 Myr of the LCC subgroup and 14 - 15 Myr of 
the UCL subgroup, respectively. These age estimates were 
revised by Sartori et al. (2003) to 16 - 20 Myr on the basis 
of the Padova isochrones (Bertelli et al. 1994) instead of 
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the Maeder (1981a, b, c) isochrones, which were used by 
de Geus et al. (1989). The Schaller et al. (1992) isochrones, 
which we used, are an upgrade of Maeder's isochrones by the 
Geneva group. Moreover we note that Sartori et al. (2003) 
have adopted for the majority of their stars the spectral 
types given in the Hipparcos catalogue, which might not be 
as reliable as the Geneva photometric data which we used. 
Given these uncertainties we conclude that our age datings 
of the LOG and UGL subgroups are consistent with Sartori 
ct al.'s result. This agrees also well with the age of Pleiades 
subgroup Bl, which was suggested to be responsible for the 
origin of the LB by BB02, but is significantly larger than 
assumed by Maiz-ApcUaniz (2001), especially for the LCC 
subgroup. Moreover, we have examined with the help of the 
Simbad data base each star of the subgroups lying not on 
the main sequence and found that practically all these stars 
are either binaries or peculiar in the sense that they are 
variable, emission hnc stars etc. (cf. the notes to the tables), 
so that their position off the main sequence in the colour- 
magnitude diagram shown in Fig. 5 can be explained in our 
interpretation by such efiiects. 

In Fig. 6 we trace back the positions of the UCL and 
LCC subgroup members over the last 30 Myr using again 
the epicycle equations (1). However, we have not used the 
individual space velocities of the stars, but adopted for each 
star the mass-weighted mean velocity of the combined sub- 
groujjs. This avoids any unphysical spread of the spatial dis- 
tribution of the stars at earlier times due to the errors of the 
space velocity components of the stars. The stellar masses 
have been determined with the mass-to-magnitude relation 

M»/M@ = 3.857 - 1.453Mv + 0.183Mv + 0.069Mv , (2) 

which we have derived from a fit to Schaller ct al.'s (1992) 
isochrone data. To the mean velocity of the stars we have 
added the solar motion {U,V,W)q = (10,5.3,7.2) km/s 
(Dehnen & Binney 1998), so that the orbits arc calculated 
in the reference frame of the local standard of rest. Since 
interstellar gas has usually only small peculiar motions, the 
local interstellar gas, and with it the LB which is indicated in 
Fig. 6 by the contour line taken from Lallement et al. (2003), 
will basically corotate with the local standard of rest around 
the Galactic centre. This means that the LB is at rest in the 
reference frame of Fig. 6. As can been seen from Fig. 6 the 
path of the association has aligned itself 15 Myr ago nearly 
parallel to the tangential ^-direction, the direction of Galac- 
tic rotation. Remnants of supernova explosions occurring 
during this period will have experienced very little shear due 
to the differential rotation of the Galaxy. The shear effect 
is described quantitatively by the term linearly proportional 
to time in the epicycle equation for Y{t) (cf. Eqs. 1). The co- 
efficient X{Q) — is the mean guiding centre radius of an 
orbit. If the spread of these radii is small, as was the case in 
the last 15 Myr, the shear effect of the Galactic differential 
rotation is minimized. In our view this might well explain 
why all supernovae occurring during that time have com- 
bined together to form the LB, while supernova remnants 
formed at earlier times have drifted away (cf. Figs. 6 and 
7). 

Next we illustrate in Fig. 7 the position of the Upper 
Centaurus Lupus and Lower Centaurus Crux associations 
relative to the LB today and at earlier times and reproduce 
the present day LB contours in meridional sections through 



Figure 6. Path of the Upper Centaurus Lupus and Lower Cen- 
taurus Crux associations over the last 30 Myr projected onto the 
Galactic plane. The look-back time is colour coded. The orbits 
are calculated backwards in the reference frame of the local stan- 
dard of rest assuming for each star the same mass-weighted mean 
velocity of the stars. The position of the Local Bubble is indicated 
by the dash-dotted contour line and is at rest in this reference 
frame. 



Figure 7. Meridional sections of the contours delineating the 
outer boundary of the Local Bubble together with the positions 
of the stars in the Upper Centaurus Lupus and Lower Centaurus 

Crux associations. The horizontal axis in the upper left panel 
points into the direction I = 300°, in the upper right towards 
I = 315° and so on. The vertical direction is always perpendicular 
to the Galactic midplane. The ages of the associations are colour 
coded as in Fig. 6. 

the bubble. From Fig. 6 we estimate the Galactic longi- 
tude in which direction we expect the associations to move. 
Choosing then the appropriate meridional section through 
the bubble from the paper by Lallement et al. (2003), we 
can determine immediately the positions of the stars in that 
longitude range relative to the LB. As can be seen from the 
upper panels of Fig. 7 the associations are today just about 
to exit the bubble. Five and 10 Myr ago they were inside. 
The bottom right panel of Fig. 7 indicates that they entered 
15 Myr ago the region occupied by the LB today. In this sce- 
nario the LB was starting to form about 15 Myr ago, which 
is consistent with the estimates of the age of the LB by 
Maiz-Apellaniz (2001) and BB02. In this context it should 
be kept in mind that although the contours determined by 
Lallement et al. (2003) are the presently best available, they 
are derived from Nal absorption line measurements, which 
allow to trace the Hi distribution under certain conditions, 
such as low temperatures (< 10'* K), line saturation for high 
column densities etc. (see Sfeir et al. 1999 for a discussion). 
In particular, the extension of the LB is unknown in di- 
rections where there are no background stars. Regardless of 
the uncertainties of the outer boundary of the LB one might 
wonder, moreover, how realistic a scenario is, in which the 
SNe explode rather close to the edge of the present day bub- 
ble. As our high resolution simulation discussed in the next 
section shows, the location of the star cluster with respect 
to the centre of the bubble is not crucial. The bubble ex- 
pands always fastest in the direction of the lowest ambient 
density and pressure. Since in the direction of the Galac- 
tic centre the Loop 1 superbubble was formed almost at the 
same time as the LB by SNe exploding within the Sco-Cen 
association, the pressure in this direction is very high. Hence 
the LB was forced to expand rather towards the anticentre 
direction and perpendicular to the plane, in agreement with 
the observations. 

3.2 Unravelling the supernovae of the Local 
Bubble 

Next wc shall derive the number and masses of stars that 
have exploded during the journey of the OB association 
through the LB. Following BB02, the number of supernovae, 
which created the LB can be estimated with the aid of the 
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(3) 



initial mass function (IMF). We have fitted an IMF of the 
form (Massey et al. 1995) 

dN _ dN 
dM^dM 

with an index F = —1.1 ± 0.1 to the data. Masses are given 
in units of solar masses. The lower end of the main sequence 
at A/v = 1 mag corresponds to AO stars with masses of 
A4i = 2.6 Al© and the upper tip at Afv = —3.7 mag to BO 
stars with masses of Mu = S.2Mq (Schaller et al. 1992), 
respectively. The total number of stars in the UCL and 
LCC associations, respectively, allow the determination of 
the normalization constants, 



N 



dN 
dM 



M'^'dM = 0.228 



dN 
dM 



(4) 



implying = 184 for the 42 UCL and l^j^^ = 118 for 

the 27 LCC stars, respectively. As we have shown in the 
previous section, OB stars entered the LB region 10 — 15 
Myr ago, setting the clock for its origin to t — 0. From a 
further fit to Schaller et al.'s isochrone data we estimate 
that the main sequence life time of such bright stars scales 
with mass as 



(5) 



with To = 1.6 • lO^yr and a = 0.932. This means that the 
masses of the most massive stars Mat in the associations 
at a lookback time of At years ago are given by 



Mat = 



To 



(6) 



implying Mio = ISAMq if At = 10 Myr or = 
26.6 Al© if At = 15 Myr depending on the entry time of 
the associations into the volume occupied by the LB to- 
day. The expected number of supernovae, i.e. the number of 
'missing' stars, is then calculated by 



NsN = 



dN 
dM 



M^'-'dM . 



(7) 



We thus obtain estimates of A'^sn =8 — 12 from the Upper 
Centaurus Lupus and Asn = 6—8 from the Lower Centaurus 
Crux associations, respectively. The estimate of 14 to 20 
supernovae, which created the LB, is in good agreement with 
the values determined by Maiz-Apellaniz (2001) and BB02. 
Extrapolating the IMF to masses beyond Mat we estimate 
that 12 to 5 SNe exploded before the associations entered 
the present LB volume. 

We have noted above that our original sample is com- 
plete in radial velocities for stars with colours {B ~ V) < 
—0.1 which corresponds to My = 0.7, if an extinction of 
Ay = 0.1 is assumed. According to Eq. (2) such stars have 
a mass of 2.95 Mq. If we remove 15 stars with colours 
(B — V) > —0.1 from the final sample (79 stars) and modify 
Eq. (4) for the cut-off at 2.95 at the low mass end, we find 
= 291 instead of 302. Thus the incomplete- 



dN 

d-M lo,UCL-t-LCC 



ness of the original sample has not introduced any significant 
bias in our sample. 

In order to assess the question, whether the estimated 
number of supernovae would suffice to excavate the LB, we 
consider the energy input by the supernova explosions into 
the interstellar gas. According to the initial mass function 
(3) there are 



dN = 



dN 
dM 



M^~^dM 



(8) 



stars in the mass range {M,M + dM) with main sequence 
life times (t, t — dr) . Thus 

dN 



dN 



dM 



M' 



dM\ ^ 

iirj''"' 



and the energy input rate is given by (cf. BB02) 



^SN = -ttSsnNsn = £sN 77— 

dt dt 



(9) 



(10) 



= £SN 



dN 
dM 



M'^-'i-l) 



d / T 

dT \T0 



dr 
dt 



where £sn denotes the energy released by a single supernova, 
£sN = 10^^ ergs. According to the way we have set up Eq. (9) 
^ is equal to 1. Equation (10) describes the trade-off of 
the increasing number of supernova progenitors and their 
increasing main-sequence life times with decreasing mass. 
Inserting the age-to-mass relation (5) into Eq. (10) leads 
then to 



iSsN = iSsNO ^7 



(11) 



with tr defined as tr = t/lO'^yr. For the constant £sno we 
find 



A £sN dN 

osm — -777 

OTo dM 



(12) 



= 3.5 • lO'' 



dN 
dM 



erg/s 



Equation IIH shows a rather weak decline of the super- 
nova energy input rate into the LB as a result of partial 
compensation between the increasing number of stars with 
decreasing mass and a corresponding increase in main se- 
quence life time. It is quite remarkable - although probably 
fortuitous - that the distribution of stellar masses during 
the star formation process is nearly anti-correlated with the 
main sequence life time of stars. 

BB02 have derived a bubble wind equation which de- 
scribes the growth of the size of the bubble with time. The 
radius of the bubble is given by 



Rhit) — RhO t ^° 

with the constant 



RhO 



475a 



(4a - 7r)(3a - 4r) 



i+r/a 



27r(2 - r)po 



(13) 



(14) 



We note in passing that Eqs. (13) and (14) are consistent 
with Eqs. (14) - (16) of BB02, except for a different value of a 
used here, and a normalization error in BB02, where Lo and 
po should be replaced by Lo = Lo/tf-^ with 5 = — (l + F/a) 
and Po = po/Rq- In Eq. (13) a constant density po of the 
ambient interstellar gas is assumed for which we adopt a 
value of Po = 2 ■ 10"^'' g/cm^. The index in Eq. (13), 
(2a — r)/5a = 0.564, lies between the index of 0.4 of the 
Sedov equation, describing supernova remnants, and the in- 
dex of 0.6 of the stellar wind/superbubble expansion law. 
For a LB age of 10 to 15 Myr Eq. (13) predicts a bubble 
radius of 78 pc to 100 pc, respectively. This in good agree- 
ment with the observed size of the LB in the Galactic disk, 
as determined by Lallement et al. (2003; cf. also Fig. 7). For 
the determination of the expected LB size we have used the 
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expected numbers of supernovae both from the LCC and 
UCL subgroups. Maiz-Apellaniz (2001) has argued that the 
LB owes its existence only to the 6 SNe stemming from the 
LCC subgroup, because stars from this subgroup came clos- 
est to the Sun in the past. We find the same when tracing 
the orbits of the stars backwards in time. However, the mem- 
bers of the UCL subgroup did enter the region occupied by 
the LB today and SNe stemming from the UCL subgroup 
have to be taken into account, in our view, in the energy 
considerations as well. The energy input of 6 SNe would ex- 
cavate a bubble with radius of only 65 pc, which is more 
difficult to reconcile with the fact that the walls of the LB 
have been blown out above and below the Galactic plane so 
that the LB has become effectively a chimney. In general, 
however, similarity solutions as applied here can only give a 
rough estimate of the LB age and size due to several severe 
restrictions. Firstly, the ambient medium has to be assumed 
to be either homogeneously distributed or to follow a power 
law distribution in density, and its pressure has to be small 
compared to the bubble pressure. Secondly, turbulent mix- 
ing and mass loading, which occur in real bubbles, are hard 
to incorporate without further assumptions (cf. Dyson et 
al. 2002). Therefore the most realistic approach to model 
existing bubbles is to perform 3D high-resolution numerical 
simulations of their formation. A first simulation of this kind 
was carried out by Breitschwerdt & Avillez (2006) which was 
based on the older and less detailed LB formation scenario 
of BB02. In the next section we present an upgrade of that 
simulation which is now based on the better understood su- 
pernova rate and the calculated paths of their progenitors 
through the LB as derived in this paper. 



3.3 High Resolution Simulations of the Local 
Bubble Evolution 

We have simulated the effects of the explosions of the stars 
formerly belonging to the UCL and LCC subgroups as 
their trajectories have crossed the LB volume towards their 
present positions. The crucial physical boundary conditions 
we have to apply to our simulations are the locations, the 
masses, and derived from this, the explosion times of the 
supernovae responsible for the origin of the LB. The latter 
can be inferred from Eqs. © and lllll to be 



dN _ 1 dN 
dr aro dM 

which can be integrated to 



VT"0/ 



AT/ N 1 dN 



(15) 



(16) 



taking the mass of the most massive star as an upper bound- 
ary. We then proceed, somewhat arbitrarily, to bin the num- 
ber of exploded stars between Mu and Mat modulo integer 
solar masses, and derive their main sequence and hence ex- 
plosion times from Eq. @. Next, the explosion locations are 
fixed by assuming that the presently "missing stars" were 
following the centres of mass of their respective subgroups. 

The 3D high-resolution simulations are based on a hy- 
drodynamical Godunov scheme (cf. Codunov & Ryabenki 
1964) supplemented by adaptive mesh refinement (AMR) 
along the lines described by Avillez & Breitschwerdt (2004) , 



Figure 8. Temperature (top panel) and pressure (bottom panel) 
distributions in the Galactic midplane 13.4 Myr after the first 
explosion in UCL occurred. The pressure is given in units of 
cm~'^ K, i.e. divided by Boltzmann's constant k. The dimensions 
and morphology of the Local Bubble are similar to the present 
observations. Loop I, to the right of the LB, is bounded by an 
X-ray illuminated shell (top panel). 



Breitschwerdt & Avillez (2006) . This entails a detailed treat- 
ment of the evolution of the interstellar gas in a volume of 
the Galaxy with a square area of 1 kpc^ and a vertical extent 
of 10 kpc on either side of the Galactic midplane based on 
the 3D SN-driven ISM model of AviUez (2000) and Avillez 
& Breitschwerdt (2004). In these calculations the ISM is dis- 
turbed by background supernova explosions at the Galactic 
rate. Initial conditions for the ambient medium were chosen 
from a data cube of a previous hydrodynamical run where 
the highest adaptive mesh refinement resolution was 1.25 
pc (Avillez & Breitschwerdt 2004, Breitschwerdt & Avillez 
2006). As a specific boundary condition we have to include 
the simultaneous evolution of the Loop I superbubble, which 
has been observed to interact with the Local Bubble accord- 
ing to ROSAT PSPC observations (Egger and Aschenbach 
1995). We therefore selected a site with enough mass to form 
all the high mass stars which are expected to explode as su- 
pernovae. Using the same IMF for Galactic OB associations 
we derived in total 81 stars with masses A4 between 7 and 
31 A4q which in our simulations compose the Sco Cen clus- 
ter; 39 massive stars with 14 ^ ^ 31A^0 have already 
gone off, generating the Loop I cavity (see Egger 1998, see 
also Avillez & Breitschwerdt 2005a). Presently the Sco Cen 
cluster, which is located at (375, 400) pc in the top panel 
of Figure |5] hosts 42 stars to explode within the next 13 
Myr. Periodic boundary conditions are applied along the 
four vertical boundary faces of our computational volume, 
while outflow boundary conditions are imposed at the top 
{z = 10 kpc) and bottom {z = —10 kpc) boundaries. The 
simulation time of this run was 30 Myr. 

Figure|H|shows the temperature (top) and pressure (bot- 
tom) distributions in the Galactic midplane 13.4 Myr after 
the explosion of the first SN, only a few thousand years af- 
ter the last UCL and LCC supernovae with masses of 8.2 
Mq have exploded. This can be seen as a red spot at {x,y) 
— (200, 300) pc. The LB is located in the region between 
100 ^ X ^ 300 pc and 250 ^ J/ ^ 550 pc, its centre being 
located at {x,y) = (200,400) pc. The shock waves of the 
last two SNe occurring within the LB are most noticeable 
in the P/k distribution by the high pressure peak shown in 
the bottom panel. To the right of the LB the shell of Loop I 
can be seen, which due to its high temperature will emit 
in soft X-rays (top panel), consistent with ROSAT PSPC 
observations. 

Another striking feature in Fig. |H] (bottom panel) are 
the coherent bubble structures within a highly disturbed 
background medium with a pressure in the range 2 5C 
log(P/fc) ^ 4 which are due to the locally enhanced SN 
rates in the vicinity of the Sun and in the Loop I region. 
The successive explosions close to the Sun heat and pres- 
surize the LB, which at first looks smooth, but develops 
internal temperature and density structures at later stages. 
About 13.4 Myr after the first explosion the LB cavity, which 
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is bounded by an outer shell will start to fragment due to 
Rayleigh-Taylor instabilities, in agreement with a linear sta- 
bility analysis carried out by Breitschwerdt et al. (2000). It 
then fills a volume roughly corresponding to the present day 
LB size. 

A more detailed analysis of these results and their ob- 
servational consequences will be the subject of forthcoming 
papers. 
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4 CONCLUSIONS AND OUTLOOK 

In contrast to previous analyses of the origin of the local 
bubble we have not merely selected presently known stellar 
subgroups and traced their kinematics back in time. Instead 
we have scrutinized ab initio a large sample volume of stars 
for stellar groups by analyzing their spatial and kinematical 
properties. From such an unbiased search among nearby B 
stars we confirm the rather robust result that besides the 
Upper Scorpius (US) subgroup the Upper Centaurus Lupus 
(UCL), and Lower Centaurus Crux (LCC) subgroups with 
ages of 20 - 30 Myr are the youngest stellar associations 
in the solar neighbourhood. Our search volume is presently 
limited to a diameter of 400 pc, because the Hipparcos par- 
allaxes are not accurate enough at distances larger than 200 
pc. Hence the analysis of a larger volume has to await the 
launch of GAIA. 

Our search strategy relied mainly on kinematical crite- 
ria, and we found many other B stars with the same kine- 
matics as the subgroups. We have followed the paths of the 
associations into the past and find that they entered the 
region of the present LB 10 to 15 Myr ago. Deriving Ovi 
column densities from a numerical simulation of the general 
ISM (Avillez & Breitschwerdt 2005b) as well as of LB and 
Loop I evolution (Breitschwerdt & Avillez 2006) in a re- 
alistic background medium, excellent agreement was found 
with Ovi absorption line data obtained with FUSE (Oegerle 
et al. 2005, Savage & Lehner 2006). According to numer- 
ical LB evolution simulations by Breitschwerdt & Avillez 
(2006), who used supernovae from the subgroup Bl of the 
Pleiades to power the LB, the Ovi data can be fitted with 
a LB age of 14.4±!J:| Myr. The age of 13.9 - 14.1 Myr es- 
timated from the present simulation is thus consistent with 
the age estimated from the slightly different simulation by 
Breitschwerdt & Avillez (2006). We therefore conclude that 
the LB must have been excavated during this time. We find 
that about 14 to 20 SNe originated from the associations 
LCC and UCL. The implied energy input into the ambient 
interstellar gas explains quantitatively the present size of the 
LB. 

The LB serves as an ideal test laboratory for superbub- 
ble models due to the wealth of observations against which 
they can be tested. Apart from the important Ovi test, we 
will also compare EUV and soft X-ray emission data with 
our models in order to derive the excitation history of ions 
in the LB and a possible deviation from coUisional ionisation 
equilibrium. 
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APPENDIX A: HIPPARCOS NUMBERS OF 
SELECTED STARS 



This paper has been typeset from a TJjX/ WFff/i file prepared 
by the author. 
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Table Al. Identified members of Upper Scorpius (ass=l), Upper Centaurus Lupus (ass=2), and Lower 
Centaurus Crux (ass=3). Positional and velocity errors are given by the e^. 
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Table Al. continued 



HIP- 


ass 


Mv 


{B - V)o 


X 


Y 


Z 


U 


V 


W 


ex 










ew 


no. 




[mag] 


[mag] 


[pc] 


[pc] 


[pc] 


[km/s] 


[km / s] 


[km / s] 














75304 
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-1.83 


-0.158 


159.8 


-78.5 


53.6 


-8.17 


-22.61 


-5.98 


22.9 


11.2 
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2.7 
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1.3 
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-0.08 
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4.6 


2.9 


76297 


2 


-3.41 


-0.217 


151.9 


-76.7 


35.8 


-8.26 


-21.78 


-8.70 


32.8 


16.5 


7.7 


8.9 


6.4 


2.9 


76371 


2 


-1.07 


-0.184 


115.1 


-63.8 


20.3 


-2.86 


-19.47 


-2.27 


11.2 


6.2 


2.0 


1.7 


1.8 


0.6 


76395 


2 


0.95 


-0.103 


99.4 


-39.3 


35.1 


-8.24 


-14.60 


-4.54 


9.3 


3.7 


3.3 


3.3 


2.0 


1.3 


76945 


2 


-0.67 


-0.145 


108.5 


-42.8 


33.6 


-3.21 


-22.09 


-3.65 


11.1 


4.4 


3.4 


1.7 


2.2 


1.0 


77286 


2 


0.23 


-0.118 


107.2 


-40.8 


32.0 


-0.87 


-19.16 


-2.89 


11.6 


4.4 


3.5 


6.7 


3.2 


2.1 
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-1.77 
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144.3 


-35.7 


59.2 


-5.12 


-19.74 


-7.55 


21.0 


5.2 


8.6 
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2.0 


77840°! 


1 


-1.31 


-0.160 


120.4 


-28.1 


49.0 


-10.45 


-15.61 


-8.52 


18.9 


4.4 


7.7 


1.5 
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1.6 
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-0.089 
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17.1 
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1.6 


2.0 


1.0 
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-0.60 


-0.140 


136.6 


-28.7 


54.1 


-12.34 


-15.49 


-10.56 


16.8 
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6.6 


3.1 


2.2 


1.5 


78265" 
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-2.85 


-0.249 


128.9 


-29.2 


48.7 


-12.31 


-15.32 


-10.47 


15.2 


3.5 


5.8 


6.8 


2.6 


2.7 


78384 


2 


-2.48 


-0.226 


1.38.4 


-53.8 


28.9 


-5.80 


-21.12 


-6.28 


16.3 


6.3 
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3.6 


2.8 


1.2 


78655 


2 


-1.14 


-0.141 


148.4 


-56.6 


29.2 


-9.36 


-22.89 


-7.30 


18.2 


6.9 


3.6 


2.8 


3.0 


1.1 


78756^5 


2 


0.83 


-0.056 


153.9 


-59.9 


28.0 


-6.36 


-22.90 


-4.40 


24.5 


9.5 


4.5 


6.9 


4.5 


1.7 


78877°* 


1 


-0.02 


-0.096 


137.5 


-23.4 


53.1 


-6.35 


-19.40 


-10.64 


17.6 


3.0 


6.8 


3.4 


2.7 


1.7 


79044 


2 


1.15 


-0.077 


120.5 


-40.4 


25.9 


-2.31 


-20.04 


-5.47 


12.5 


4.2 


2.7 


6.9 


3.0 


1.7 


79404 


1 


-1.32 


-0.197 


134.4 


-28.3 


41.6 


-5.18 


-15.43 


-6.86 


15.2 


3.2 


4.7 


6.9 


2.4 


2.3 


81914 


2 


0.15 


-0.119 


141.9 


-44.0 


8.1 


-6.73 


-14.66 


-3.91 


18.0 


5.6 


1.0 


1.8 


2.0 


0.6 


82545 


2 


-2.47 


-0.223 


153.6 


-37.7 


10.7 


-3.26 


-19.78 


-3.70 


20.9 


5.1 


1.5 


0.9 


2.7 


0.8 


84970 


1 


-2.95 


-0.223 


171.6 


1.4 


19.7 


-1.37 


-20.10 


-5.05 


20.5 


0.2 


2.3 


3.6 


2.3 


1.0 



'^multiple, -^emm.l. /variable, ^variable, ^rotnl. variable, ^Be, ^emm.l. /binary, '^Be/neb. emm., 
^variablc//3Ccp, ^^binary, ^^ccl. binary, ^'^variable/ZJCep, ^^spec. binary, ^^double, 
variable / /3Cep /double? , ellips. variable / double? 



Table A2. HIP-nos. of the remaining selected stars 



2484 


2505 


5566 


7588 


7943 


8886 


10602 


10944 


11249 


12692 


15338 


15404 


15444 


15627 


15770 


15988 


16147 


16210 


16244 


16470 


16611 


16803 


17499 


17531 


17563 


17573 


18033 


18190 


18213 


18216 


18788 


19860 


20042 


20063 


20171 


20186 


20884 


21192 


21281 


22109 


23607 


23767 


24244 


25813 


26248 


26487 


26623 


26640 


29426 


30069 


30122 


30675 


31278 


31362 


31685 


32677 


32912 


33015 


33579 


34045 


35054 


35785 


36188 


37304 


38455 


38863 


39138 


39360 


39906 


40581 


42177 


42637 


43105 


43394 


43878 


43937 


45080 


45418 


45941 


46283 


47119 


47391 


47452 


51437 


51576 


52370 


52419 


52502 


52701 


52736 


54767 


55597 


58484 


60000 


60718 


60855 


61789 


62786 


65474 


66821 


67301 


68269 


69389 


71353 


76243 


76669 


78493 


79653 


82673 


82902 


83895 


85391 


85792 


86414 


89482 


89908 


90200 


90422 


91235 


91729 


92614 


92855 


93104 


93187 


93231 


93299 


93805 


95400 


95951 


96052 


96417 


96468 


97376 


98412 


98754 


100751 


101017 


101421 


101475 


101716 


101746 


101868 


103089 


103532 


103616 


104105 


105148 


105282 


106604 


107462 


107664 


107930 


108022 


109139 


112781 


115990 


116805 
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